Stationary flow fields prediction of variable physical domain based on proper orthogonal decomposition and kriging surrogate model

In this paper a new flow field prediction method which is independent of the governing equations, is developed to predict stationary flow fields of variable physical domain. Predicted flow fields come from linear superposition of selected basis modes generated by proper orthogonal decomposition(POD). Instead of traditional projection methods, kriging surrogate model is used to calculate the superposition coefficients through building approximate function relationships between profile geometry parameters of physical domain and these coefficients. In this context,the problem which troubles the traditional POD-projection method due to viscosity and compressibility has been avoided in the whole process. Moreover, there are no constraints for the inner product form, so two forms of simple ones are applied to improving computational efficiency and cope with variable physical domain problem. An iterative algorithm is developed to determine how many basis modes ranking front should be used in the prediction. Testing results prove the feasibility of this new method for subsonic flow field, but also prove that it is not proper for transonic flow field because of the poor predicted shock waves.     

应用傅里叶时间谱方法求解二维跨音速流动问题解的精度研究

为了从时间离散精度的角度评估傅里叶时间谱方法的实际计算效率,对用该方法求解二维跨音速流动问题的解进行了精度研究。通过求解Euler方程模拟俯仰振动NACA0012翼型的周期性非定常流动。结果表明:傅里叶时间谱方法可以使用很大的时间间隔较为准确地模拟有激波存在的周期性跨音速流动;对于此类流动问题,相比二阶向后差分公式法,傅里叶时间谱方法在预测翼型表面压力系数积分量特别是升力系数的时间变化规律方面具有明显的计算效率优势;用于模拟翼型表面压力系数自身的时间变化规律时,傅里叶时间谱方法具有不低于二阶向后差分公式法的计算效率。     

Shock protection of penetrator-based instrumentation via a sublimation approach

It is often necessary for space-borne instrumentation to cope with substantial levels of shock acceleration both in the initial launch phase, as well as during entry, descent and landing in the case of planetary exploration. Current plans for a new generation of penetrator-based space missions will subject the associated on-board instrumentation to far greater levels of shock, and ways must therefore be found to either ruggedize or else protect any sensitive components during the impact phase. In this paper, we present an innovative method of shock protection that is suited for use in a number of planetary environments, based upon the temporary encapsulation of said components within a waxy solid which may then be sublimated to return the instrument back to its normal operation. We have tested this method experimentally using micromachined silicon suspensions under applied shock loads of up to 15,000g, and found that these were able to survive without incurring damage. Furthermore, quality factor measurements undertaken on these suspensions indicate that their mechanical performance remains unaffected by the encapsulation and subsequent sublimation process.     

A self-similar solution of a curved shock wave and its time-dependent force variation for a starting flat plate airfoil in supersonic flow

The problem of aeroelasticity and maneuvering of command surface and gust wing interaction involves a starting flow period which can be seen as the flow of an airfoil attaining suddenly an angle of attack. In the linear or nonlinear case, compressive Mach or shock waves are generated on the windward side and expansive Mach or rarefaction waves are generated on the leeward side. On each side, these waves are composed of an oblique steady state wave, a vertically-moving one-dimensional unsteady wave, and a secondary wave resulting from the interaction between the steady and unsteady ones. An analytical solution in the secondary wave has been obtained by Heaslet and Lomax in the linear case, and this linear solution has been borrowed to give an approximate solution by Bai and Wu for the nonlinear case. The structure of the secondary shock wave and the appearance of various force stages are two issues not yet considered in previous studies and has been studied in the present paper. A self-similar solution is obtained for the secondary shock wave, and the reason to have an initial force plateau as observed numerically is identified. Moreover, six theoretical characteristic time scales for pressure load variation are determined which explain the slope changes of the time-dependent force curve.     

挂滴燃烧中的微汽泡行为和液滴跳动

利用高速显微摄像技术,观察研究了RP-3微尺度挂滴燃烧过程中,微汽泡的核化、生长、聚并和溢出的行为过程及液滴跳动现象.结果得出:①液滴中微汽泡产生的位置;②在蒸发初期时由微汽泡引起液滴体积的膨胀;③在燃烧中由汽泡行为导致液滴的跳动;④液滴燃烧的稳定(不爆裂)现象.分析了表面张力在挂滴燃烧过程中的作用,确定了挂滴稳定燃烧的临界条件下的临界特征体积在2.5~5.0μL.      

航天器火工冲击技术研究进展

系统总结了国内外航天器火工冲击环境预示、地面试验和减冲击技术三个方面的研究进展,并分析了国内在以上三个方面与国外航天强国的差距。在此基础上,从我国航天工程实际需求出发,分析指出了今后航天器火工冲击领域应重点开展的研究方向。     

Numerical investigation of the impact of asymmetric fuel injection on shock train characteristics

Numerical simulations are carried out to investigate the impact of asymmetric fuel injection on shock train characteristics using the commercial-code FLUENT. The asymmetry of fuel injection is examined by changing the fuel flow rates of the upper and lower wall fuel injectors. The numerical approach solves the two-dimensional Reynolds-averaged Navier–Stokes (RANS) equations, supplemented with a k-ω model of turbulence. As a result, different ways of fuel injections will always lead to shock train transitions, with the variations of shock train structure, strength and leading edge position. For symmetric fuel injection, the flowfield of the isolator is quite asymmetric with the boundary layer of the upper wall side developing much stronger than that of the lower wall, which is due to the heterogeneity of the incoming flow. Regarding to asymmetric fuel injection with more of lower wall side, though the pressures in the combustor are nearly the same, the first shock of the shock train converts between ‘Distinct symmetric X type shock’ and ‘Obscure and weaker asymmetric shock’ and the shock train leading edge moves upstream with the increase of the asymmetry level. With regard to asymmetric fuel injection with more of upper wall side, ‘incomplete asymmetric X type shock’ occurs and the shock train structures keep nearly the same with low level of fuel injection asymmetry. Unexpected results like unstart will happen when increasing the level of fuel injection asymmetry. And the isolator will come back to normal state by decreasing the differential of upper and lower wall sides fuel injections.     

Solar cycle variation of “killer” electrons at geosynchronous orbit and electron flux correlation with the solar wind parameters and ULF waves intensity

To construct models for hazard prediction from radiation belt particles to satellite electronics, one should know temporal behavior of the particle fluxes. We analyzed 11-year variation in relativistic electron flux (E>2 MeV) at geosynchronous orbit using measurements made by GOES satellites during the 23rd sunspot cycle. As it is believed that electron flux enhancements are connected with the high-speed solar wind streams and ULF or/and VLF activity in the magnetosphere, we studied also solar cycle changes in rank order cross-correlation of the outer radiation belt electron flux with the solar wind speed and both interplanetary and on-ground wave intensity. Data from magnetometers and plasma sensors onboard the spacecraft ACE and WIND, as well as magnetic measurements at two mid-latitude diametrically opposite INTERMAGNET observatories were used. Results obtained show that average value of relativistic electron flux at the decay and minimum phases of solar activity is one order higher than the flux during maximum sunspot activity. Of all solar wind parameters, only solar wind speed variation has significant correlation with changes in relativistic electron flux, taking the lead over the latter by 2 days. Variations in ULF amplitude advance changes in electron flux by 3 days. Results of the above study may be of interest for model makers developing forecast algorithms.     

Variations of the solar acoustic high-degree mode frequencies over solar cycle 23

Using full-disk observations obtained with the Michelson Doppler Imager (MDI) on board the Solar and Heliospheric Observatory (SOHO) spacecraft, we present variations of the solar acoustic mode frequencies caused by the solar activity cycle. High-degree (100 < ? < 900) solar acoustic modes were analyzed using global helioseismology analysis techniques over most of solar cycle 23. We followed the methodology described in details in [Korzennik, S.G., Rabello-Soares, M.C., Schou, J. On the determination of Michelson Doppler Imager high-degree mode frequencies. ApJ 602, 481–515, 2004] to infer unbiased estimates of high-degree mode parameters ([see also Rabello-Soares, M.C., Korzennik, S.G., Schou, J. High-degree mode frequencies: changes with solar cycle. ESA SP-624, 2006]). We have removed most of the known instrumental and observational effects that affect specifically high-degree modes. We show that the high-degree changes are in good agreement with the medium-degree results, except for years when the instrument was highly defocused. We analyzed and discuss the effect of defocusing on high-degree estimation. Our results for high-degree modes confirm that the frequency shift scaled by the relative mode inertia is a function of frequency and it is independent of degree.     

Fe/O Ratios in Interplanetary Shock Accelerated Particles

It is widely accepted that diffusive shock acceleration is an important process in the heliosphere, in particular in producing the energetic particles associated with interplanetary shocks driven by coronal mass ejections. In its simplest formulation shock acceleration is expected to accelerate ions with higher mass to charge ratios less efficiently than those with lower mass to charge. Thus it is anticipated that the Fe/O ratio in shock-accelerated ion populations will decrease with increasing energy above some energy. We examine the circumstances of five interplanetary shocks that have been reported to have associated populations in which Fe/O increases with increasing energy. In each event, the situation is complex, with particle contributions from other sources in addition to the shock. Furthermore, we show that the Fe/O ratio in shock-accelerated ions can decrease even when the shock is traveling through an Fe-rich ambient ion population. Thus, although shock acceleration of an Fe-rich suprathermal population has been proposed to explain large Fe-rich solar particle events, we find no support for this proposal in these observations.